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Abstract 



The temperature dependence of quark and gluon condensates in QCD 
as precursor of the chiral and deconfining phase transition is modelled 
O l' with a conformally extended non-linear a-model including broken chiral 

CD . and scale invariance. The model is further enlarged by including (free) 

■ heavier hadrons. Within this frame we then study the interplay of QCD 

scale breaking effects and heavier hadrons in chiral symmetry restoration. 

X 

H ' 1 Introduction 

Quantum Chromodynamics (QCD) with light quarks possesses, at the classical 
level, the remarkable feature of chiral and conformal symmetry. At low en- 
ergy, these symmetries are broken by the non-perturbative vacuum structure 
which is signalled by the appearance of nonvanishing quark and gluon conden- 
sates, (0|<7q|0) and (0|Cr" G° |0), respectively Q. As a result of spontaneous 
breakdown of chiral symmetry there arise light Goldstone bosons (pions) which 
are the relevant hadronic degrees of freedom of low-energy QCD. On the other 
hand, at higher temperature the non-perturbative vacuum structure partially 
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disappears which manifests itself in a phase transition from hadronic matter 
with confined quarks and gluons to an interacting plasma of deconfined quarks 
and gluons. In particular, it is important to investigate the temperature be- 
haviour of the quark condensate in order to determine the critical temperature 
T c at which the chiral symmetry becomes restored. This question has been ex- 
tensively studied in lattice QCD Q from first principles. A natural question 
in the context of particular models is what influence the existence of heavier 
(non-Goldstone) mesons a, p, ui, ... and baryons N, A, ... has on the actual 
chiral symmetry restoration temperature T c . Finally, taking into account also 
the temperature-dependent gluon condensate, it is interesting to consider the 
interplay of temperature effects on the condensates and other thermodynamic 
quantities. 

For the study of non-perturbative condensates we use in this paper the 
effective Lagrangian approach. Effective meson Lagrangians provide a compact 
and extremely useful method to summarize low-energy theorems of QCD ||. 
They incorporate the broken global chiral and scale symmetries of QCD. In 
particular, to mimic the QCD scale anomaly 0] 



KZ D ) = (^g^) , (i) 



one introduces a scalar chiral singlet dilaton-glueball field \ with an interaction 
potential V(x) (|] so that, via Noether's theorem MM cx x 4 and (0) is satified 
for xo — (0|x|0) 7^ 0. As proposed in || fjj the kinetic and symmetry breaking 
mass terms of standard chiral meson Lagrangians have then to be multiplied 
by suitable powers of x/Xo m order to reproduce the scaling behaviour of the 
analogous terms of the QCD Langrangian. These ideas have incited a number 
of investigations of the chiral and scaling behaviour of low-energy Lagrangians 
at finite temperature §, g, [h], [n], |§ [l| . 

Here, we will start with a non-linear (scaled) 0(4) cr-model which is isomor- 
phic to a SU (2) x SU (2) chiral Lagrangian. The model contains one scalar (a) 
and three pseudoscalar if fields, constrained by the scaled condition a 2 + if 2 — 

(xo) wnere ft i s tne P^ on d eca y constant. In the path integral evaluation 
of the thermodynamic potential this constraint can be rewritten by introduc- 
ing a Lagrange multiplier field X(x). This allows to perform the remaining 
Gaussian integration over the pion fields exactly and to apply a saddle point 
approximation to the A integral. Finally, to model the temperature effects of 
heavier hadrons, we add the free Lagrangian of these particles. As discussed in 
[pL~4| (without studying the interplay with the gluon condensate) their influence 
would lead to a lowering of T c by about 10 percent. Therefore, a generalization 
of the model which includes the effect of the gluon condensate on the low-lying 
hadronic states (with mass, say, < 2 GeV) is worthwile. To couple these non- 
Goldstone particles to the remaining degrees of freedom in our model, we have 
to adopt a simple recipe defining the scaling properties of their masses with 
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X- In the spirit of Ref. |14| , we use a simple-minded additive quark model to 
determine the chiral-symmetric and non-symmetric pieces of hadron masses. 
This leads to a definite scaling prescription for hadron masses Mh — * Mh[x/Xo\- 
Notice that the finite temperature glueball dynamics might influence the chiral 
symmetry restoration, too |ll| . In order to clarify this point and to compare 
the thermal effect of glueballs with that of the other heavy hadrons, we include 
quantum fluctuations of the glueball field. The effective masses M x and M e ff 
of glueballs and pions are then obtained in a self-consistent way from saddle 
point equations derived for the thermodynamical potential (free energy) of the 
model. 

The structure of the paper is as follows. In sect. 2 we introduce the effective 
low-energy meson Lagrangian with broken chiral and scale invariance discussing 
in some detail the links between chiral- and scale-symmetry breaking on the 
basis of the low-energy theorem for the gluon condensate in the presence of 
non-vanishing quark masses jl5| . Moreover, the scaling behaviour of the masses, 
which is different for Goldstone and non-Goldstone particles, is set up. In 
sect. 3 we calculate the free energy density of the extended non-linear cr-model 
including pions, heavier hadrons and glueballs and determine the temperature 
dependent quark and gluon condensates ((qq))\r a (T) (including the effect 
of heavier hadrons) and ((G 2 ))|t °c X 4 from the saddle point equations for the 
a and \ fields, respectively. Finally, sects. 4 and 5 contain the discussion of 
numerical results and the conclusions. 



2 The effective hadron Lagrangian 

2.1 Symmetry breaking and scale anomaly 

In order to model the approximate chiral and scale invariance of QCD we con- 
sider an extended 0(4) non-linear cr-model containing one scalar a- and three 
pseudoscalar if -fields as well as a scalar isoscalar dilaton (glueball) field \- The 
corresponding Lagrangian is given in Euclidean notation by 

L(a, If, x) = \{d,af + \{d^f + \{d, X ) 2 + V( X ) + V SB (a, X ) , (2) 
where the fields satisfy the rescaled chiral constraint 

with /tt = 93 MeV being the pion decay constant, and where xo = (0|x|0) has 
been introduced for dimensional reasons. This is the actual vacuum expectation 
value of the field x which will slightly depend on the amount of explicit breaking 
of chiral symmetry. 
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In order to take into account the gluon contribution to the QCD scale 
anomaly we have included the dilaton interaction potential S 

V( X ) = K x 4 ( log ( 2L )- \ ) (4) 



Xq) 4 

which takes its minimum value at X — Xq- Finally, Vsb(&, x) is the scaled chiral 
symmetry breaking term of scaling dimension 3 

\ 2 

X 



V SB (x,v) = -cv[ ^ 1 (5) 

to be added to V(x)- It should be noted that due to the constraint (||) the field 
a = cr(x, If*) is not an independent degree of freedom. For instance, omitting 
at tree level the pion field (giving rise to the pion loop contributions), we have 
a(x, ) = fn^- The total potential 

Vtot(x) = V(x) + VsB(x>*(x>~ft)) ( 6 ) 
has then its minimum shifted to 

Xo=X Q + 2-J-3 c U + 0(c 2 ). (7) 

^ Xq 

The parameters K , \Q an d c will be specified later by considering the vacuum 
energy density (bag constant), the glueball mass and the pion mass. 
From eqs.(||||!]) we get the trace anomaly 

- -K X 4 - cfJX-X , (8) 



Xa 



and the vacuum energy density 

Cvac = \{%i f )=V tot ( X0 ) 

= -\k x 4 o-\cU (9) 

= -\Kx\ - cf n + 0(c 2 ), 

where we have used eq. (^|) in the last line. 

It is instructive to compare this with the QCD expression for Nj = 2 

= y~ G%G% +m(uu + dd) , (10) 
z 9 
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where j3 denotes the Gell-Mann-Low function, 

and where m — m u — md are the light current quark masses. [] 

Taking into account the low-energy theorem for the gluon condensate in the 
presence of non-vanishing quark masses fl5|| 

( Yf^ G lv > = ( Y~ g G ^ G % )|m=o + 3m(uu + dd) (12) 

which relates the shift of the gluon condensate compared to the massless case 
to the quark condensate, one obtains for the vacuum energy density the form 

= \ (®%Z D ) = { ^G^G^) lm=0 + m (uu + dd) . (13) 

The comparison with eqs. (^) and ( |l3| ) suggests the following identifications: 

( G%G%, )\ m=0 = -K( x 4 )|c=o = - K X 4 q (14) 

and 

m ( uu + dd ) = - c A = - Ml fl . (15) 



Eq. ([L5|) is just the well-known Gell-Mann-Oakes-Renner relation. We have 
already used that, discarding pion loops, (a) = oq = f„ due to the constraint 
(||). Moreover, by expanding in eq. (|^) the 7? dependent field a 



a=u f f-wsr (16) 

to leading order in if 2 one finds that the pion mass, given by M n = \J c/ f n , 
scales as 

— » M n ,[-^- . (17) 
V Xo 



2.2 The different scaling behaviour of masses 

One of the main issues of the present work is to study the interplay of the quark 
and gluon condensates at finite temperature (near the chiral phase transition) 
and their relation to the hadron spectrum. In particular, we are interested in the 
question how the light u and d quark condensate and the gluon condensate will 
be affected by the inclusion both of additional massive hadrons h = p,uj,N, ... 

x We do not discuss effects due to the anomalous dimension of the gq-operator. 
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(represented by bosonic or fermionic fields iph) and a (lowest) glueball state. 
Note that for the inclusion of non-Goldstone particles we need to know also the 
scaling properties of the corresponding mass term. Using the relations |l6[ 

(H(k)\Q^\H(k))=2M* , (18) 

(N(p)\e^\N(p)) = M N «(p) u(p) (19) 

(for bosons h — H and fermions h — N) and taking eq. (|l0| ) into account, we 
see that the masses Mh contain chiral-symmetric and non-symmetric pieces. A 
more detailed investigation of this issue requires the use of low-energy theorems 
of QCD (l^, [l7|] which is outside the scope of this paper. Here we shall follow a 
procedure which determines the hadron masses from a simple-minded additive 
quark model as follows ^ 

M h = N% m u + N£ m d (20) 

where N% is the number of valence quarks (or antiquarks) of type a = u,u,d 
and d with masses m a inside hadron h. rh a denotes the total quark mass of 
quark type a which decomposes as 

m a = m a + m at d yn (21) 

into current and dynamically generated mass. Let, for illustration, the dy- 
namical part of the quark masses vn a ^ yn be represented by a scale-invariant 
interaction term of the form 

Lint = -gaqq~-gf 7T — qq = - m dy n — qq (22) 

Xo Xo 



where use has been made of the approximate expression of eq. ( |16[ ) (omitting 
the pion field) and of the Goldberger-Treiman relation mdyn = 9 fir- From eq. 
( p2|) we deduce the scaling behaviour of the dynamical quark mass 

Wldyn * mdyn • (23) 

Xo 

In the case of Goldstone bosons (pions) one has the mass formula 

Ml oc m m d yn (24) 

which formally reproduces the scaling behaviour of eq. (0). For the non- 
Goldstone particles, on the contrary, we obtain from eqs. ( pp| , pl| , p3| ) the 
scaling prescription 

M h M h \x/Xo] (25) 

with some well-defined functions Mh of the ratio x/xo (Mh should take for 
X = Xo the known value Mho)- In order to include non-Goldstone hadrons 

2 For similarity and differences compare with Ref. [11.4]. 
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we shall simply use a free Lagrangian L h (ifh, Mh[x/xo}) for each particle type. 
This corresponds to the assumption that interaction terms with pions lead to 
higher order derivative terms in the effective Lagrangian of the pion sector. Such 
terms will be neglected in thermodynamical considerations at low temperatures 
studied here. The same holds for interactions between non-Goldstone particles 
because of their low thermal density below the critical temperature. 

3 Chiral and conformal thermodynamics 

3.1 The free energy density of the extended non— linear 
a— model 

In the following we shall study the thermodynamical properties of the system of 
pions, glueballs and free non-Goldstone hadrons described by the Lagrangian 

L = L(a, t, x) + L ' l (^ M h\x/Xo]) ■ (26) 

h 

The partition function of the Lagrangian (|2^) reads 



' DX Dy,{<T,l?,x,Vh 

J e— z oo J 




(27) 



+ 2 ( fl fX) +V tot (a, X ) 

+ YL h (cp h ,M h \x/xo})jj 

where Du(...) denotes the integration measure of the physical fields. j3 = 1/Tis 
the inverse temperature, and we have rewritten the 5-function constraint 8{a 2 + 
if 2 — ,f 2 (-^) 2 ) following from eq.(^) in terms of an integral over the auxiliary 
field X(x). Periodic (antiperiodic) boundary conditions for boson (fermion) fields 
as well as the appropriate form of the particles' Lagrangians L h ((ph, Mh\x/xo\) 
are understood. In the following we shall evaluate the path integral in eq. ( P7| ) 
in the saddle point approximation for A and a (putting X(x) = X = const, 
<j{x) = a — const and dropping the corresponding integrations) but shall keep 
the quantum fluctuations of the glueball field x( x ) = X + x( x ) 0- We expand 

3 Note that the saddle point approximation in A corresponds to the leading order in a 1/N 
expansion of the 0(N) cr-model MM . 
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the integrand up to terms 0(x 4 )- Performing first the Gaussian integration over 
the fields if and <fh we obtain (x = (r, x)) 



J' J Dx exp(-J*dr J d 3 x($[A,a,x] + ~(^x) 

. d<& i _ 2 a 2 $ i 9 3 $ i _ 4 a 4 $\ . 



with 




+ ^t ot ((T,x) 



2 ( l0 S(- 92 + 2A ))(^) (M) 



3 

]T 9h {\og(-d 2 + M h [x/xo] 2 )) {x , x 

esons h 

E 9h (tr D irac log (9 + M h [x/xo]) 



baryons h 



Here, gh are statistical degeneracy factors (g p = g n = 4, g u = 3,g p = 9 etc.). 
The first logarithmic term in eq.(p9|) leads just to the partition function of the 
free relativistic Bose gas of pions with an effective (temperature dependent) pion 
mass 

M e 2 // = 2 A. (30) 

In order to get rid of the term linear in x in eq. (Pq) let us choose x to satisfy 
the following saddle point equation 

§ = -2 A ft * + ^ ± |:£^r lo g (-M[ 2) [^]) = 0. (31) 
OX Xo »X ox^ Xo 

The last sum refers to bosonic (fermionic) hadrons and the arguments ... of 
the functional logarithms stand for their respective wave operators. Moreover, 
to simplify further calculations, we shall take into account the interactions of 
massive hadrons with the glueball field x onr y U P to first order, discarding 
(on the level of first order perturbation theory) also the term 0(x 3 )- In this 
approximation nucleons do not contribute to the coefficient of the 0(x 2 ) term 
due to their Yukawa-type coupling to the glueball field x> anc l the massive 
hadron contribution in the coefficient of the 0(x 4 ) term can also be discarded. 
On the other hand, the interaction term -gx 4 (9 4 Vtot(o',x)/d x 4 ) leads, in the 
mean field approach, to a dynamical contribution to the glueball mass. To see 
this, it is convenient to rewrite the 0(x 4 ) term by introducing an auxiliary field 
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p giving rise to a Yukawa-type coupling p\ 2 • Indeed, applying a usual Gaussian 
integration we have 



exp 



(-jf*/*^**) 







DpeXP {~J dT j d ^\-^ t p2+l 2 P * 2 )) (32) 
where the fourth derivative of the potential is given by 



V t { :> (x) ^2AK log (J-j + 44 K, (33) 

and the factor 3 in the l.h.s. of eq. ( p^ ) accounts for the three possible pairings 
to express x 4 as a square of % 2 . Inserting the identity ( |32| ) into eq. (|2^), we 
read off an expression for the effective glueball mass ^] 

w (-**(-' 

dx \ VXo, 



M 1(P) = sI-^It) + V tot (a,x) 



n ]\,f2 

Y,9h^ n(M ft 2 [x/xo]) + P 



= l2Kx 2 \o g {^\ + 4Kx 2 - ^1 _ £^ (34) 
\Xq/ Xo Xo 



o 7i/r2 

^ Xo 



ft 

Here, II(M 2 ) is the expression for a closed meson loop given by the sum of T ^ 
and T = contributions 

n(Af 2 ) = n T (M 2 ) + n (Af 2 ) , 
d 3 P 



with £(p) = (p 2 + M 2 )a . Inserting eq.(p2[) into eq.(p8[) we can perform the 
integration over the glueball field x- Thus, there remains a functional integral 
over the p field, 



J Dp expl-fdr J (-f + l(lo g (- d 2 + Ml(p))) ix M(36) 



9 V} 

, z v tot 

4 We use here ^ 2 r~/^„l ~ /i™ 2 (— u,„\l ^ (— u,..\ 2 



M llx/xo] « 4m^ n (x/xo) ~ A^ (x/xo) neglecting small terms O(m), 
0(m 2 ) depending on the current mass. 
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which will be evaluated at the saddle point solving the gap equation 

Po = n(M^W). (37) 

Finally, we have 

Z = exp(-mF eff (a,\,x)) (38) 
where V3 is the 3D volume and the free energy density F e f f is given by 

F e ff{a,X,x) = V tot ((T,x) + A ^' j2_ ^(^) j + ( continued ) 

+ F*(M* ff ) + FX(M*(p )) (39) 
- 1 Vg} (n(Af x 2 ( Po ))) 2 + ]T F h (M h [x/xo}) ■ 

h 

Here F v , F x and F h are the expressions for the thermal determinants including 
the zero point energy [] 

F*,x,h = F *' x > h + F^ x ' h (40) 

with 



effJ 



(27T) 
d 3 p 



Fml) = r/^tog^-e-VS^/T), (42) 
J#(M fc ) = ± 5 ,T|^log(l T e-V^/T 

« -9H^-e-^' T . (43) 

The upper (lower) signs in eq.(^) correspond to bosons (fermions) with the 
masses Mh substituted by M/Jx/Xo] depending on % Finally, Fq, Fq and Fq 
denote the zero temperature contributions to the free energy density of the pion, 
glueball and hadron fields, 

A - d 3 p 1 
(2tt) 2 2 



A * d 3 P 1 

(2tt) 2 2 



A, 



d 3 p 1 



F h (M h ) = ± 5 ,/ ^- V / P 2 + M 2 . (46) 



5 Note thatn(M 2 ) = 2g^ y . 
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Also here the masses Mh are understood depending on x/xo- Note that, due 
to the necessary momentum cut-off, a naive regularization (with fixed cut-off) 
would violate the scaling properties of our effective theory which is constructed 
to substitute QCD. In order to keep the wanted scaling behaviour, we have 
introduced a rescaled cut-off A x = following the analogous argumentation 
in the case of the Nambu-Jona-Lasinio model jnj, 19 1. Moreover, to avoid 
additional parameters, we specify the cut-off in the spirit of effective chiral 
theories by A = 4 n f n = 1.2 GeV for all fields. 

Following Ref.Q, we now adopt a suitable subtraction procedure which 
guarantees the following two renormalization conditions at T = 0, 

M^Ao) = 2 A = Ml = -f , (47) 

( v ) = o-o = f w . (48) 

Here Ao and o~o are the known values of A and a which should extremize F e ff 
at T = 0. In addition, here we have to fulfill a third renormalization condition 

( X) —Xo ■ (49) 

For this purpose let us define a subtracted expression of the free energy den- 
sity F su b(a, A, x) which is obtained from F e ff(a, \,x) ( |39"| ) by replacing in the 
determinantal terms F n — > F^ ub with 



= FZ{Mi ff )+FZ{Mi }f 



dFl_ 

^effJ | T=0 



and using analogous subtraction prescriptions for F x , F h and for the expressions 
n = 2^jy f\ . The subscript T = means that the derivatives have to be 
evaluated at the physical zero temperature masses, i.e. at M e ff = M x = 139 
MeV for the pions, for the non-Goldstone hadrons at Mh = Mho a s known from 
the Rosenfeld table, and at M x = M g b being the mass of the lowest glueball 
state. Note that in eqs.(|34|) and ([37]) II(Af 2 ) has to be replaced by the subtracted 
expressions H su b(M 2 ). 

Then, the thermal averages, abbreviated as <j{T) = (( a ))t, A(T) = (( A ))x 
and x(F) = (( X ))tj are obtained from solving the following saddle point 
equations 

$F su fy 3F su b dF su b , . 

-do~= ' — = ° ' -W = Q ' (51) 



e Strictly speaking, one has also to introduce subtraction terms — (A x — A) ^ (with 

zero temperature masses Mo). Due to the weak dependence of A x on T these terms turn out 
to be numerically very small. We neglect them putting A x — * A. 
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taking the derivatives at a = cr(T), A = A(T) and x = x{T)- Note that in eq.(pT{) 
we have replaced the "classical" value x by a function x{T) which satisfies a 
saddle point equation including now the contribution of the glueball loop. Q 
The first two equations take a simple form 

2 ^^>-((l?) S + sl%) < 52 » 

and 

[ ) U Xo ^ II X(T) 2 d{Ml ff ) + x{T?d{MD) ■ [b6) 

No massive (non-Goldstone) hadrons except the glueballs contribute here be- 
cause their mass does not depend on A. Clearly, eq.(|5^ ) is just the thermal 
average of the constraint (|J) evaluated in the saddle point approximation. No- 
tice also that the glueball loop contributes here with a different sign than the 
pion loop. 

It is instructive to consider for a moment the chiral limit c = 0. Then cq.([52"|) 
admits the two possibilities, 

case (i) A(T) = , cr(T) ^ (54) 

and 

case («) A(T) ^ , a(T) = . (55) 

Clearly, case («) describes spontaneous breakdown of chiral symmetry with a 
nonvanishing order parameter er(T). One expects that this holds for tempera- 
tures T < T c , where T c is the critical temperature, with vanishing A(T) reflecting 
the existence of a massless Goldstone pion. On the contrary, case (ii) should 
hold for T > T c where chiral symmetry is restored and pions become massive. 
The situation should be qualitatively similar also for c ^ 0. 

Inserting the solutions <r(T), A(T) and x(T) into F su b yields the actual free 
energy density as a function of temperature 

F(T) = F sub (a(T),X(T), X (T)) . (56) 

By definition, the subtracted quantity F(T = 0) is just the vacuum energy 
density (||), associated to the gluon and quark condensates through eq.(|l~3|). 
Indeed, using X (T = 0) = X o, A(T = 0) = A = <t(T = 0) = /„ and 

c = f n M%, we obtain 

F(T = 0) = V tot (U, X o) 

= -\ K xt--Mfl (57) 



4 ™ 4 
-Kx A q -Ml fl + 0{Mt) 



7 For a justification of this procedure in terms of quantum (loop) corrections to the effective 
action we refer to Refs. 



t tins pro( 
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We will now determine the parameters K and xo (Xq) m cq- (|57|) taking the bag 
constant B and some lowest glueball mass M g b at T = as two reference scales 



B = - W = -K X t + -Ml Si (58) 



4 4 



and 



m2 d 2 VtM,x) ft 
M gb = ^=2 — 59 

M 2 f 2 
Xo 

For the bag constant we choose Bi = 240 MeV, which is compatible with 
the gluon condensate at T = extracted from QCD sum rules For the 
glueball mass we take M x — 1.6 GeV, a value motivated by lattice results 
Experimental searches favour candidates in the mass region of 1.5. ..1.8 GeV 

3.2 Quark and gluon condensates at T ^ 

Let us next investigate the influence of the glueball field and of the heavier 
(non-Goldstone) nonstrange hadrons on the chiral condensate of u and d quarks. 
Primarily, the chiral condensate is obtained as the logarithmic derivative of the 
partition function with respect to the current quark mass m, 

II- \\ T dlogZ dF(T) 

((qq))T = -T?—F. — = —^ — ■ 60 

V am dm 
With eq.(|T5|) we have 

d ( uu + dd )n d ( uu + dd 



(61) 
(62) 



dm d c fl dMl ' 

such that 

(( qq )) T _ 1 F(T) 

(qq)o U d c 

where use has been made of the relation 
d F(T = 0) 

« = ~U ■ 63 

o c 

Calculating (( qq ))t we have to take into account also the c-dependence of 
the glueball mass M 2 (|34|) and of the heavy hadron masses in eq.(^5|) via the 
current quark masses ( |20| ) and (^l|). From eqs.(|6^) and (^3|) we obtain 



{{qq)) T <?(T) f X (T) 



qq )o A V xo 



2 



A , (64) 
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where A contains contributions from the glueball and the other massive hadron 
states 



A = — 



2 <j(T) d Ff ub 1 ^ d Ml d F* ub (M h ) 



U xl d(M$ f v V 9c 5(M < 
with Mh = Mh[x(T)/xo]- Here we have used 
d (Af 2 -pa) 2 a 



f_ Ft r Ft (M?\ ' 1 ' 



and 



d c xl 



dM MT)/Xo] _ 2Mh[x{T)/xo] dM MT)/xo] 



(66) 



o c a c 

-W + Nfj ._ 2 ^^. , (67) 

( uu + dd ) \ Xo / 

where and iV^ are the numbers of light valence u, d quarks (or antiquarks) 
in the hadron h. We postpone the discussion of the last term A in eq.(|64|) and 
consider the first one at first. Using eq.(p33|) and discarding there the glueball 
loop term in the bracket as well we find 



«TC»T _ ( X(T) V /j _ * 



* (rn d-*-*- i m 



Xo 7 V 2/2 8/4 



with 



Xl d F: ub 



p AO sub (f\<X\ 

H - X {Tf d A |A=A(T)- ^ 
For the special case c = 0, when M e //(T) = Af^ = (putting A x = A), we get 

Notice that if we set x{T)/xo ^ 1 we recover the structure of results obtained in 
chiral perturbation theory fllij with inclusion of one- and two-loop contributions 
[| In the case c ^ our loop expressions contain a thermal pion mass M e ff(T) = 
2 A(T) which has to be determined numerically in a selfconsistent way from 
eqs.@||,||). 

8 The two-loop contribution oc _R 2 //* in eq. (^) differs from the corresponding term in the 
SU(2) X SU(2) matrix formulation numerically. Ibis deviation is due to next— leading order 
terms in the 1/N expansion not considered here 
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Let us consider next at which critical temperature T c the order parameter 
(( QQ ))t vanishes such that chiral symmetry is restored. From eq.(|68|) we obtain, 
assuming x(T c ) ^ 0, that T c is determined by 

R{Tc) = fl • (71) 

Moreover, in the simplified case of massless quarks (c = 0, m = 0), we can use 
cq.(|70|) to derive the estimate (\o = X<j) 

T c = 2 U = 186 MeV . (72) 

Xo Xo 

The factor x(^c)/xo can be expressed by the gluon condensate. Indeed, taking 
the thermal average of eqs. d^JTo|) we obtain for m — 0, c = 

-K (( X 4 »t = (( ^- )) T , (73) 

or 



(( G 2 )) T (( X 4 »t ^ /x(T) 



G 2 )o Xo V Xo 



4 



(74) 



where G 2 = G" G" and the approximation (( x 4 ))t ~ (( X ))t = x(7") 4 has 
been made. From these considerations we expect 

t - = ( { \^-1\ )t ) 4 186 Mev ■ ( 7s ) 



G 2 ) 

As we shall find below, (( G 2 ))t/( G 2 )o < 1 in the region where (( qq ))t ^ 0. 
The gluon condensate depends only very weakly on T in this region. There- 
fore, the estimate 2 f v for the chiral transition temperature remains practically 
unchanged by it. Let us next consider the positive glueball term in eq . (|53|) . Al- 
though it is small in comparision with the pion loop, its contribution is expected 
to be non-negligible in the region T ~ 2 where the first two terms in eq. ( j53[) 
approximately cancel, shifting thus T c upwards. On the other hand, as was ar- 
gued in Ref. (l4) (without consideration of gluon condensation and glueballs), 
the inclusion of massive hadron states alone leads to a negative contribution to 
the quark condensate and gives rise to a lowering of T c by about 10 percent. In 
order to estimate the strength of this effect in our model we have to consider 
the additional terms in eq.(p4) omitted until now. Notice first that the arising 
glueball contribution in eq.(|65|) is - in addition to the exponential suppression 
with the high mass - suppressed by the small factor er(T) compared to the con- 
tribution of other heavier hadrons and can thus be safely neglected. Clearly, 
the remaining negative hadron contributions will have a tendency of lowering 
T c and thus counteract its increase originating from the positive glueball term 
in eq.(|53|). The numerical investigation of the interplay of both effects will be 
done in Section 4. 
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Finally, the temperature dependence of the gluon condensate ( |74| ) is numer- 
ically determined from the third saddle point equation dF su b/9X\ x (T) = m 
eq.(|5l|) together with eqs. (|52|) and (|53|). 



4 Discussion of numerical results 
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Figure 1: Temperature dependence of the partial pressures Pi(T) for pions, for 
heavier mesons, for baryons and the contribution from the dilaton potential 

V( X )- 

In this section we present the results of the solution of the saddle point equa- 
tions. In Fig. [l] we show for sake of illustration, the temperature dependence of 
the thermal pressure(s) p%(T) of pions, all heavy mesons and baryons as well as 
the tree (potential) term, for the chiral limit m — > 0. Since the a model taken 
here describes only Nf = 2 light flavours, we have included only non-strange 
heavier hadrons. 

In Fig. the thermal average c(T) is presented together with the quark 
condensate as function of temperature in the chiral limit. One finds that the 
gluon condensate (included to mimic the breaking of scale invariance in QCD) 
has only a negligible influence on the behaviour of cr(T) and on the quark con- 
densate. The vanishing of <j{T) at some T c is dominated by the intrinsic chiral 
dynamics of pions modified by an increase of T c by about 10 MeV (6 percent) 
due to glueballs, but is not related to the small reduction of x{T) with rising 
temperature. This conclusion is in accord with Ref. |ll[] and does not confirm 
the corresponding conjecture of Ref. 0. We see that heavier non-Goldstone 
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Figure 2: Temperature dependence of the thermal average er(T) of our model 
(compared with the cr-model for pions only) and the quark condensate (((?<z))t 
for the chiral limit using a bag constant i? 1 / 4 = 240 MeV. All quantities are 
normalized to their values at T = 0. 

hadrons indeed influence the quark condensate in the region of T < 186 MeV 
due to the current quark mass dependence of the hadron masses. They yield 
an effective lowering by about 10 MeV of the increased value of T c obtained 
for the cr-model with pions and glueballs leading to a restoration of the value 
T c w 2/jr = 186 MeV obtained for pions alone. These results are in qualitative 
agreement with the conclusions of Ref. E3| where pions and heavier hadrons 
have been included and contributions to the thermodynamical potential have 
been calculated up to three loop level for pions (while no interplay with the 
gluonic degrees of freedom has been considered). 

Fig. U exhibits the same quantities for the case of finite quark masses m = 
m u = rrid « 7 MeV (explicit chiral symmetry breaking) which are adjusted to 
give the pion a mass = 139 MeV. Without QCD scale breaking effects taken 
into account a is now shifted upward throughout all temperatures and does not 
play anymore the role of an order parameter. With the coupling to the gluon 
condensate and, through the glueball mass, to the glueball loop, a even bends 
upward for temperature T > 220 MeV. This can be interpreted as a kind of 
stabilization of a by the gluonic degrees of freedom. A similar behaviour has 
been found in Ref. [Q. 

If we neglect the smaller glueball loop contributions we see that the tem- 
perature dependence of a(T) is directly related to the ratio x(T)/M e ff(T) (sec 
eq. (p2h. Finally, including the negative contributions of the other heavier 
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Figure 3: Same as Fig. 2 for a finite pion mass M„ = 139 MeV (to = 7 MeV) 
using the same bag constant. 
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Figure 4: Temperature dependence of er(T), the effective pion mass M e ff(T), 
the gluon condensate ((GG))t and the glueball mass M X (T) for the chiral limit 
using the same bag constant as before. The effective pion mass is normalized 
to = 139 MeV, all other quantities to their respective values at T = 0. 

hadrons, the total quark condensate is again strongly decreased and vanishes at 
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Figure 5: Same as Fig. 4 for M n = 139 MeV (m = 7 MeV). 

T c ps 220 MeV. In Fig. [| we show the temperature dependence of the effective 
pion mass M e ff(T), the gluon condensate ((G 2 ))|t oc \ A an d the effective glue- 
ball mass M X (T) in the case of the chiral limit. All quantities are normalized 
to their respective values at T = 0, with the only exception of Me//, which is 
always normalized as M e ///M w , whether c — > or not. Notice that for T > T c 
there arises a non-vanishing pion mass from the cr-model, whose teperature 
behaviour is now influenced by the gluon condensate. The plot ends when A 
hits the constraint A > 0. 

Fig. |^ presents the same quantities as Fig. 4 in the case of finite quark 
and pion mass. Clearly, for temperatures T > 180 MeV the above discussion 
and the figures must be taken with caution with respect to real QCD. From the 
point of view of the effective hadronic Lagrangian one has to recall that in this 
temperature region higher order derivative terms have to be included into the 
action and higher loop contributions are expected to contribute. Therefore our 
results are at best qualitative ones in the temperature region T > T c . 

5 Conclusions 

In this paper we tried to give a qualitative description of various important as- 
pects of temperature-dependent QCD using an effective Lagrangian approach. 
Our analysis is based on a conformally extended non-linear cr-model which de- 
scribes light pions coupled to the gluon condensate and, eventually, to glueball 
thermal loop effects. Most importantly, this model embodies explicit and dy- 
namical breaking of chiral and conformal symmetry by construction. Heavier 
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non-Goldstone mesons and baryons have further been added in order to model 
the additional dependence of the condensates via heavy hadron masses. 

The primary motivation was to investigate in this framework the mutual 
influence of quark and gluon condensates. However, if the parameters are fixed 
using information from zero temperature like the vacuum energy density and 
the glueball mass, we find that the gluon condensate is rather temperature inde- 
pendent in this physical frame up to temperature T m 200 MeV. At somewhat 
higher temperature, in the interval up to T — 250 MeV, it decreases to 0.75 of 
the zero temperature value. This effect is accompanied by a (roughly inverse) in- 
crease of the glueball mass. If the bag constant is taken smaller [B = (140MeV) 4 
instead of B = (240MeV) 4 ) this "melting" and the stronger rise of the glueball 
mass with temperature set in at slightly lower temperature. 

In the chiral limit, the temperature dependence of the a field, the quark 
condensate etc. is almost independent of the gluon condensate, but slightly 
changed by the glueball loop contributions. It is strongly influenced by the 
explicit chiral symmetry breaking through the quark mass term. Let us first 
summarize the case of the chiral limit (m u = md = = 0). Discarding 
loop effects of glueballs and heavy hadrons, the critical temperature of chiral 
symmetry restoration turns out to be T c = 2 / OT x ^ ( f ~ 186 MeV, in view of 
X(T c )/xo ~ !• The inclusion of the glueball loop in a(T) then increases this 
value of T c by about 10 MeV. Next, heavy hadrons and dynamical glueballs 
contribute with different signs to the extra term A (cf. eq. (|65|)) in the quark 
condensate. Due to the large glueball mass and the small value of the factor a(T) 
near T c the respective glueball term is here numerically too small to counteract 
the decrease of the quark condensate (and of T c ) caused by heavy hadrons. 
In accord with Ref. O], we find a 6 percent (10 MeV) downward shift of T c 
which just compensates the corresponding increase arising from u{T) in our 
model. In particular, our results confirm the conclusion of Ref. jll| that in the 
chiral limit the vanishing of the quark condensate is, above all, dictated by the 
internal pion dynamics. This behaviour is only slightly corrected by glueballs 
and the heavier non-Goldstone mesons and baryons and has nothing to do with 
a decrease (not to mention vanishing !) of the gluon condensate. Thus, the 
conjecture of Ref. (?]] is not supported by our analysis. Note also that at T > T c 
the effective pion mass becomes non-zero, grows to a maximum and decreases 
to zero again. This behaviour is related to the beginning decrease of the gluon 
condensate in the temperature interval 200 MeV < T c < 250 MeV. However, 
beyond the temperature of chiral symmetry restoration the predictions of the 
cr-model (present here in its simplest version) have to be taken with caution 
and are at best qualitative ones. 

In the more realistic case of soft explicit breaking of chiral symmetry by finite 
quark and pion masses, the curve of the pion loop contributions to the quark 
condensate (represented by the saddle point solution <r(T)) is shifted upwards 
and does not clearly define a phase transition. (This is the case even for the 
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pure a model.) Discarding the small glueball loop contributions, this behaviour 
can be directly related to the behaviour of the ratio x(T)/M e ff(T) as can be 
seen from eq. (|52|). 

Nevertheless, including also non-Goldstone hadrons, their contribution is 
sufficient to make the quark condensate vanish, however at somewhat higher 
temperature. This happens at T c » 220 MeV. Again, we have to emphasize 
that this temperature region is beyond the region of applicability of the simple 
non-linear cr-model used here as a corner stone of our extended model. From 
chiral perturbation theory []l4| |24| it is known that higher order derivative terms 
in the Lagrangian and higher loops cannot be neglected at these temperatures. 

Summarizing, our model exhibits in the chiral limit, for temperatures T < 
200 MeV, competing effects from glueballs (increase of T c ) and heavier hadrons 
(decrease of T c ) which, because of cancellation, have however only a minor 
influence on the chiral symmetry restoration. The gluon condensate shows also 
very weak temperature dependence in this temperature range. On the other 
hand, in the more realistic case of finite quark and pion masses, it is just the 
quark mass dependence of the heavier hadron masses which leads to a significant 
lowering of the quark condensate and, finally, to its vanishing. 

Moreover, noticeable variations of the gluon condensate down to 75 percent 
of the zero temperature value are only seen in a temperature region 200MeV < 
T c < 250 MeV where the simple non-linear cr-model is not realistic anymore. 

Modelling the implicitely temperature dependent heavy hadron masses by 
another dependence on the order parameters might lead to stronger and yet 
more interesting to analyse effects for the detailed temperature dependence of 
chiral symmetry restoration. Such models need dynamical input from hybrid 
models including quark degrees of freedom simultaneously with the Goldstone 
fields. These mechanisms will be the subject of future studies. 
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